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The  contract  was  initiated  under  Project  No.  73^2,  "Fundamental 
Research  on  Macromolocular  Materials  and  Lubrication  Phenomena" , 

Task  No.  73^203#  "Fundamental  Principles  Determining  the  Behavior 
of  Mac romole cule s" .  It  was  administered  under  the  direction  of  the 
Air  Force  Materials  Laboratory,  Research  and  Technology  Division, 
with  Dr.  Ivan  J.  Ooldfarb  as  project  engineer. 

This  report  covers  work  conducted  from  June  1963  to  May  1964. 

The  assistance  of  Dr.  Maurice  King  in  the  mathematical  analysis 
is  gratefully  acknowledged. 


The  parameter  TQ  in  the  empirical  equation 
log  f)/ d  -  log  kq  ♦  By(T-T0) 

is  being  investigated  as  an  indicator  of  glass  formation  in  simple  liquids 
for  which  viscosity  9  and  density  d  have  been  published  at  various  temper¬ 
atures  T.  The  parameter  B*  was  found  to  equal  Tp  +  330*  for  polystyrene 
fractions,  but  this  relation  does  not  hold  for  the  n -alkane  series.  Some 
correlations  have  also  been  discovered  between  the  above  constants  and 
those  in  the  Antoine  equation  for  vapor  pressure  of  the  same  liquid. 

It  was  found  that  the  Gibbs -DlMarzlo  theory  correctly  predicts  the 
values  of  T0  of  the  n-alkanes  with  a  stiffness  energy  of  490  cal  ./mole. 

The  predicted  T0  for  polyethylene  is  12$*%,  which  corresponds  to  a  well- 
known  transition.  A  detailed  study  of  sterlc  factors  in  polymers  is 
expected  to  lead  to  methods  of  predicting  the  stiffness  energy  and  hence  T0. 

A  free  volwe  theory  has  been  developed  which  accounts  satisfactorily 
for  the  change  in  viscosity  with  external  pressure,  and  which  separates 
the  phenomenological  effects  of  temperature  and  pressure  from  the  structural 
and  molecular  factors  involved  in  flow  behavior.  The  adjustable  constants 
have  not  yet  been  correlated  with  chemical  structure,  however. 

This  technical  documentary  report  has  been  reviewed  and  is  approved. 
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SYMBOLS 


t)  viscosity 

d  density 

pre •exponential  constant  in  equation  (1) 

B^  energy  constant  in  equation  (1) 

T0(i))  lover  temperature  limit  for  liquid  state  from  equation  (1) 
V  molar  volume 

h  Planck's  constant 

N  Avogadro's  number 

A?*  activation  free  energy  in  Eyrlng  equation 

R  gas  constant 

E  energy  of  vaporization 

I  structural  constant  in  Souders'  relation 

M  molecular  weight 

p  vapor  pressure 

Ap  pre -exponential  constant  in  equation  (4) 

Bp  energy  constant  in  equation  (4) 

Tq(p)  lover  temperature  limit  for  liquid  state  from  equation  (4) 
x  number  of  chain  atoms  in  polymer 

Ss  entropy  of  disorder  per  chain  atom 

Z  coordination  number 

7/  valence  of  chain  atom 

v  free  volume  fraction 

g  exp  (-£/RT) 

6  energy  difference  between  rotational  isomers 
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occupied  volune 
free  volune 

constant  In  equation  (9) 
constant  (0.2058)  in  equation  (9) 
internal  pressure  factor  in  equation  (9) 


1.  Introduction 


Molecular  theories  of  the  liquid  state  sre  at  present  Incapable 
of  quantitatively  predicting  the  effects  of  structure  of  the  molecules 
on  the  formation  of  glasses,  and  reliance  must  be  olaced  on  empirical 
relations.  The  glass  temperature  has  been  defined^)  as  that  temper¬ 
ature  at  which  the  relaxation  tine  becomes  long  compared  with  the 
duration  of  the  experiment.  KauanannV1)  also  pointed  out  that  there 
Is  a  close  similarity  between  the  relaxation  processes  which  give  rise 
to  the  change  In  liquid  structure  with  temperature  or  pressure,  and 
the  molecular  processes  In  viscous  flow.  It  was  therefore  of  interest 
to  find  that  the  best  empirical  equation  for  the  effect  of  temperature 
on  viscosity  contains  a  parameter  with  the  phenomenological  character¬ 
istics  of  a  glass  transition.  The  equation  is 


log  t) 


log  ♦ 


B»j/  (T* 


T0), 


(1) 


apparently  first  proposed  independently  by  Fulcher 
Hesse'3]#  although  it  is  often  referred  to  as  the  ' 


(2) 


,  and  Tammann  and 
Vogel  equation.  The 

T0  has  the  nhysicai  significance  of  a  lower  limit  to  the 
liquid  state.  Miller'1*'  has  pointed  out  that  the  WLF'5J  equation, 
which  is  widely  used  to  represent  the  effect  of  temperature  on  visco¬ 
elastic  behavior  of  polymers.  Is  equivalent  to  equation  (1) .  Hence  TQ 
can  be  calculated  from  the  reported  WLF  constants  for  polymers .  Values 
have  also  been  reported  for  a  large  number  of  low  molecular  weight 

liquids '3  #0,7,0)  # 


In  spite  of  this  large  number  of  determinations  o**  TQ,  no  correlation 
with  molecular  structure  has  yet  been  discovered.  Part  or  the  problem 
lies  in  the  fact  that  the  three  parameters  in  equation  (1)  are  highly 
correlated,  vhlc.  complicates  the  problem  of  calculating  the  best  values 
from  the  experimental  data.  Some  of  the  reported  T0  values  may  be 
unreliable,  due  to  the  use  of  measurements  made  over  a  small  temperature 
range  far  above  To'  and  improper  weighting  of  the  data.  Ideally,  the 
correlation  with  structure  should  be  derived  theoretically,  and  the 
adjustable  parameters  determined  by  direct  comparison  with  experimental 
data. 


Eyring  et  al.^9)  arrived  at  the  following  expression  for  the 
viscosity  of  liquids,  using  the  theory  of  absolute  reaction  rates: 

ln(V^)  -  ln(hN)  ♦  AfVrT  (2) 

where  V  ■  molar  volume,  N  *  Avogadro's  number,  h  «  Planck's  constant 
and  Af*  is  the  activation  free  energy.  It  was  also  proposed  that 
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AF*»  E/2.45,  where  E  is  the  energy  of  vaporization.  Equation  (2) 
is  not  a  good  representation  of  the  effect  of  temperature  on  viscosity, 
since  it  fails  to  predict  the  rapid  increase  in  viscosity  as  the 
glass  temperature  is  approached.  Furthermore,  some  trial  calcu¬ 
lations  of  the  ratio  E/AF*  indicated  that  this  is  highly  variable. 
Nevertheless,  the  theory  Justifies  the  comparison  of  viscosity 
behavior  with  vaporization  data,  and  such  correlations  will  be  dis¬ 
cussed  below. 

Another  correlation  which  has  enjoyed  some  success  is  the 
Souders'^OJ  relation: 

log  (log  rj  )  «  Id/M  -  2.9  (3) 

where  I  is  a  viscosity-constitutional  constant,  calculated  from 
structural  values  which  Souders  has  tabulated.  This  equation  was 
tested  against  Doolittle* s(H)  data  for  trldecane  and  heptadecane, 
and  the  API*12'  data  for  pentane.  The  plots  were  indeed  linear, 
except  at  high  densities,  but  did  not  show  a  common  intercept  at 
d  «  0.  Souders  also  claimed  the  same  slope  for  all  linear  alkanes, 
but  this  was  not  confirmed  either.  The  equation  is  not  very  useful 
as  it  stands,  but  it  indicates  one  approach  to  a  method  of  correlation 
with  structure. 

Equation  (1)  is  of  the  same  form  as  the  Antoine  equation  for 
vapor  pressure: 

log  P  ■  log  Ap-Bp/Oivr,,).  (It) 

Comparison  of  the  constants  in  equations  (1)  and  (4)  for  the  liquids 
for  which  both  have  been  reported  revealed  that  TQ  is  roughly  the  same 
for  both  vapor  pressure  and  viscosity,  but  B_  is  several  times  larger 
than  B^.  This  i6  reasonable,  since  more  energy  is  required  to  remove 
a  molecule  from  a  liquid  than  to  simply  move  it  to  a  new  location 
within  the  volume  of  liquid.  It  would  be  enlightening  to  have  a 
comparison  of  these  parameters  for  a  homologous  series,  to  see  if 
they  iollow  the  same  trends  with  structure.  Furthermore,  since  the 
energy  of  vaporization  E  can  be  calculated  from  Bp  and  T  (p)  (assuming 
ideal  vapor  behavior  and  neglecting  the  liquid  yolune  relative  to 
the  vapor),  correlations  of  E  with  structure 1*3)  could  be  translated 
into  correlations  for  viscosity  parameters. 

Preliminary  work  with  reported  viscosity  data  indicated  that 
kinematic  viscosities  obey  equation  (1)  to  ever  better  precision  than 
absolute  viscosities.  Statistical  tests  are  required  for  confirm¬ 
ation  of  this  observation. 


? 
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2.  Effect  of  Temperature  on  Viscosity 


a.  n-Paraffina 

The  effect  of  molecular  structure  on  physical  properties  can 
be  established  by  searching  for  a  function  which  accurately  represents 
the  data  for  a  homologous  series  of  materials.  For  example,  Small(13) 
found  the  parameter  (EV)1/*  to  be  linear  in  number  of  carbon  atoms 
for  several  homologous  series  of  liquids.  The  coefficient  equals  the 
group  constant  for  a  -CH2-  group,  and  the  intercept  is  that  for  the 
functioned,  or  end  group.  Since  the  parameters  in  equation  (1)  had 
not  been  reported  for  any  homologous  series,  it  was  decided  to 
determine  them  for  the  n-paraffin  series,  for  which  extensive  data 
are  ava liable » 12 ' . 

The  Mathematical  Analysis  Group  was  requested  to  calculate  ti.* 
best  values  of  log  A,.,  TQ  and  for  each  compound  for  which  data  was 
available.  At  the  same  time,  the  use  of  kinematic  instead  of  absolute 
viscosity  in  equation  (1)  was  investigated.  The  relation  between  the 
Antoine  constants  and  those  in  equation  (1)  was  also  examined. 

Table  1  gives  the  values  of  the  parameters  for  both  absolute  and 
kinematic  viscosity  derived  from  the  data  of  reference  12  by  the  use 
of  a  non-linear  least  squares  method,  programmed  on  the  Datatron 
computer.  Unfortunately,  these  data  had  been  pre-smoothed  according 
to  sorre  unknown  scheme,  and  therefore  no  statistical  tests  can  be 
made  cf  the  precision  of  the  estimates.  The  residual  mean  square  is 
substantially  lower  for  kinematic  viscosities  of  all  the  compounds 
except  butane .  As  will  be  evident  later,  the  data  for  butane  are  out 
of  11 n*  in  other  respects  as  well. 

Table  <;  gives  values  of  log  A^,  and  T0  for  both  kinematic  and 
absolute  viscosities,  calculated  from  Doolittle's'11)  data.  95^ 
Confidence  limits  for  the  parameters  are  also  included.  Again  the 
kinematic  viscosity  is  represented  more  closely  by  equation  (1). 

For  a  comparison  of  Antoine  constants  with  those  in  equation  (1), 
vapor  pressure  data  reported  by  Willingham,  et  al.'12*)  were  treated 
by  the  same  non-linear  )p«st  oquares  method.  The  results  are  shown 
,in  Table  3  along  with  the  T0  values  derived  from  the  kinematic 
viscosities  of  reference  12.  The  95^  confidence  intervals  for  the 
estimates  of  the  parameters  of  the  Antoine  equation  given  in  Table  3 
are  notably  small  since  they  do  not  include  things  such  as  an 
opera tor-to-opera tor  or  laboratory- to-iaboratory  component  of  error. 
Hence,  estimates  of  these  constants  from  other  equally  precise  sets 
of  data  might  very  well  seem  not  to  be  in  agreement  because  of  these 
other  variance  components.  The  comparison  of  T0(p)  with  T0 (*7 )  suffers 
accordingly,  and  it  would  be  necessary  to  obtain  very  precise  viscosity 
data  in  order  to  decide  whether  the  two  values  are  actually  identical. 
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ESTIMATES  OF  PARAMETERS  HI  EQUATION  (1)  FOR  THE 
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ESTIMATES  OF  PARAMETERS  Hi  EQUATION  (1)  WITH  95^  CONFIDENCE  LIMITS 
_ BASED  OH  RAW  DATA  ON  (SELECTED)  HYDROCARBONS* 
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It  was  expected  from  the  Eyring  hole  theory (9)  that  the  para¬ 
meters  B  and  6u  would  be  closely  related,  since  they  are  related  to 
the  activation  Energies  for  vaporization  and  flow,  respectively. 
Figure  1  shows  the  correlation  plot  for  (kinematic  viscosity) 
from  Table  1  vs.  B  listed  in  Table  20k  for  reference  12. 


The  trend  in  TQ  with  increasing  number  of  C  atoms  shown  in 
Tables  1  and  2  is  somewhat  irregular,  making  it  difficult  to 
extrapolate  to  the  limiting  value.  Since  this  is  the  quantity  we 
are  primarily  interested  in  correlating  with  structure,  an  effort 
was  made  to  locate  viscosity  data  on  higher  polymers.  Tung'1^) 
reported  zero  shear  melt  viscosities  as  a  function  of  temperature 
for  branched  and  unbranched  polyethylene  fractions,  which  should 
be  suitable  for  calculation  of  T0.  However,  the  points  were  too 
scattered  for  a  reliable  estimate. 


Aggarwal,  et  al.(^)  reported  viscosities  at  different  temper¬ 
atures  for  low  density  polyethylene  fractions,  and  these  were  fitted 
to  equation  (1).  Two  of  the  sets  of  data  were  too  scattered  to  be 
useful,  and  TQ  for  the  other  three  ranged  from  ^  0  to  >220*K. 

Marker,  et  al.(^)  gave  melt  viscosities  for  a  whole  polymer  of 
low  density  polyethylene  with  Mn  =  33»330 -  The  data  fit  equation  (1) 
with  Tc  «  33*K  and  =  2,642*K.  These  values  do  not  seem  to  fit  in 
with  the  results  for  the  n-paraffin  series.  It  was  concluded  that 
the  available  viscosity  data  on  polyethylene  are  not  sufficiently 
accurate  for  our  purposes,  possibly  due  to  shear  degradation  or 
oxidative  degradation. 


While  this  work  was  in  progress,  the  paper  by  Karapet* yants 
and  Yan'lci)  appeared,  giving  the  constants  in  equation  (1)  for  35 
hydrocarbons.  However,  these  workers  found  it  necessary  to  use  a 
separate  set  of  constants  for  the  upper  (>0.d5  cp)  and  lower 
(<LQ.b5)  ranges  of  values,  making  a  total  of  six  adjustable  para¬ 
meters  for  each  liquid.  Furthermore,  many  of  their  T0  values  are 
negative,  in  contrast  to  those  in  Tables  1  and  2.  all  of  which 
(except  for  butane)  are  positive.  Other  workers^;^)  have  also 
reported  negative  values  of  TQ.  This  seems  to  result  from  the  use 
of  viscosity  data  at  temperatures  far  above  T0,  and  possibly  from 
improper  weighting  of  data. 

b.  Polystyrene 

The  viscosity  data  of  Fox  and  Flory^^)  on  polystyrene  fractions 
car.  r>e  used  to  calculate  log  A^,  B*  and  T~  as  a  function  of  degree 
of  polymerization,  P.  The  cata  naa  already  beer,  fitted  to  an 
equation  of  the  same  form  as  (1)  by  Williams'*^',  and  as  Miller^ ) 
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Figure  1.  From  Smoothed  Kinematic  Viscosities  vs.  B  From  Ref.  12. 
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has  pointed  out,  the  constants  In  equation  (1)  can  be  obtained  by 
a  simple  transformation.  The  resulting  values  of  Ba  and  TQ  are 
shown  plotted  against  P  in  Figure  2.  Both  B«  and  T^  approach  limit 
ing  value 6  above  P  «  200,  and  drop  off  rapidly  below  that.  The 
strong  correlation  between  these  constants  is  evident  from  Figure  2 
There  is  an  almost  constant  difference  of  350*  between  and  T0. 

3 .  Correlations  with  Structure 

A  simple  correlation  between  energy  of  vaporization  E  and 
chemical  structure  was  proposed  by  Small'^),  based  on  the  observa¬ 
tion  that  the  function  (EV)V2  iS  linear  in  number  of  repeating 
unite  in  a  homologous  series.  Small  published  a  table  of  group 
values  from  which  E  can  be  calculated  for  any  liquid,  provided  the 
molar  volume  V  is  known.  This  relation  wa6  checked  against  the 
selected  values  of  E  and  V  for  the  n-alkanes  published  by  API 
Project  44”  2),  an(j  VQE  indeed  found  to  be  obeyed  remarkably  well, 
although  the  slope  and  intercept  were  slightly  different  from  the 
values  recommended  by  Small.  Several  drawbacks  to  this  scheme 
should  be  pointed  out,  however.  First  is  the  requirement  of  having 
a  measured  value  for  the  density  of  the  liquid.  Secondly,  the 
correlation  is  known  to  be  valid  only  at  25*C .  and  has  not  been 
tested  at  other  temperatures.  Third,  the  most  serious  objection 
is  that  the  molar  volume  is  a  strong  function  of  x,  the  number  of 
carbon  atoms,  so  that  errors  in  prediction  of  E  tend  to  be  observed 
on  a  plot  of  (EV)1/2  VB.  x. 

By  analogy,  it  was  thought  that  (BpV)1/2  should  be  linear  in 
x,  but  this  is  not  the  case.  It  has  been  found,  however,  that 

(Bpv)°-77  *  i,609  x  ^  155  (5) 

represents  the  n-alkane  data  very  well. 

Another  function  was  found  which  gives  the  effect  of  structure 
(in  the  n-alkane  series)  on  the  Antoine  constants  with  good 
precision,  and  does  not  require  knowledge  of  the  density: 

(To/Bp)2,3  «  1.12  •  10"\c  (6) 

This  result  Illustrates  the  fact  that  the  parameters  in  equation 
(4)  are  highly  correlated,  and  hence  it  should  be  possible  to 
develop  a  relation  with  only  one  or  two  adjustable  parameters. 

The  published  data  on  n-alkanes  have  been  submitted  to  the 
Mathematical  Analysis  Group  in  order  to  ascertain  the  best  values 
of  the  constants  in  equations  (5)  and  (6)  and  their  confidence 
intervals . 
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The  purpose  of  this  line  of  approach  is  to  find  a  sieans  of 
relating  intermolecular  forces  to  structure,  and  use  this  in  the 
prediction  of  glass  temperatures.  The  product  B_V,  for  example, 
is  roughly  proportional  to  "a"  in  the  van  der  Waals  equation  and 
hence  can  be  used  as  a  measure  of  cohesive  energy.  It  would  be 
highly  desirable,  however,  to  liberate  this  scheme  from  the  require¬ 
ment  that  the  molar  volume  be  known.  Future  work  will  be  directed 
toward  the  estimation  of  Bp  for  other  classes  of  compounds,  and  its 
relation  to  structural  features  rather  than  the  molar  volume . 

The  two  correlations  with  structure  discussed  above  ^equations 
(5)  and  (6 ]7  relate  to  the  vapor  pressure.  Efforts  to  extend  theoe 
to  the  viscosity  parameters  in  equation  (1)  have  not  led  to  any 
useful  result.  However,  a  rough  correlation  has  been  found  between 
Tq(i^)  and  Bp.  The  ratio  Bp/T0(»))  equals  approximately  19  for 
straight-chain  aliphatic  liquids  above  butane,  but  varies  over  a 
wide  range  for  other  types  of  compounds,  as  shown  in  Table  4. 


Table  4 


CORKEIATION  OF  Bp  WITH  ^(9) 


Liquid 

Liquid 

Bp/Tof*}) 

Methane 

8.68 

Ethylene  bromide 

18.99 

Ethane 

25.44 

Acetic  acid 

20.24 

Propane 

28.63 

Butyric  acid 

20.46 

Butane 

-73.33 

Carbon  Tetrachloride 

17.02 

Pentane 

19-53 

Cyclohexane 

24.46 

Hexane 

19-89 

Benzene 

15.89 

Heptane 

17.79 

Toluene 

107.52 

Octane 

19.47 

Chlorobenzene 

44.86 

Nonane 

17.60 

Nitrobenzene 

21.58 

Decane 

18.91 

Aniline 

9-58 

Undecane 

18.41 

Phenol 

9-63 

Dodecane 

19.15 

1 ,6  -D  ime  thylnaphthalene 

14.05 

Tridecane 

19.17 

Te trade cane 

19.72 

Pentadecane 

19.50 

Hexadecane 

19.97 

Hep tade cane 

19-41 

Octadecone 

19.92 

Nonadecane 

19.31 

Eicosane 

19.31* 
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The  values  in  Table  4  were  calculated  fro®  the  data  given  by 
the  American  Petroleun  Institute for  Bp  of  the  n -alkanes,  the 
T0(h)  values  (kinematic  viscosity)  for  these  llouida  listed  in 
Table  1,  the  values  of  Bp  given  by  Dreisbach^J  for  the  liquids 
in  colunn  two,  and  the  Tq(m)  values  for  these  reported  by  Hoffmann 
and  Rother(o).  These  are  practically  all  the  liquids  for  which 
both  Bp  and  Tp(^)  have  been  reported.  Further  investigation  along 
this  line  will  require  evaluation  of  raw  data  on  viscosities  or 
vapor  pressures. 

This  correlation  is  potentially  valuable,  since  it  relates 
T0(rj)  to  a  parameter  (Bp)  which  can  in  turn  be  correlated  with 
structure  through  equation  (5). 


The  empirical  correlation  discussed  above  are  probably  not 
sufficiently  reliable  for  extrapolation  to  infinite  molecular 
weight.  Free  volume  theories  are  also  of  little  use  for  this 
purpose,  but  the  statistical  mechanical  theory  of  the  glass 
transition  proposed  by  Gibbs  and  DiMarzio'“J  is  on  sufficiently 
firm  grounds  so  that  some  degree  of  confidence  can  be  placed  in 
its  predictions.  The  equations  given  in  reference  (22)  yield  the 
following  expression  for  the  entropy  of  disorder  per  atom,  £s, 
when  a  chain  of  x  atoms  of  valence  Z/  is  placed  on  a  lattice  of 
coordination  number  Z: 


(7) 


where  v  is  the  fractional  free  volune,  g  ■  exp(-d/HT)  and  d  is  the 
energy  difference  between  rotational  Isomers.  The  point  at  which 
As  decreases  to  zero  is  assumed  here  to  equal  T0(i}) .  Application 
of  this  theory  to  glass  temperatures  of  polymers  yields  reasonable 
values  of  6 ,  although  there  are  no  independently  measured  values 
for  comparison.  In  order  to  evaluate  the  theory  more  critically, 
the  values  of  T0(ty)  given  in  Table  1  for  the  n-alkanes  were  compared 
with  the  predictions  of  equation  (7) .  For  these  liquids,  d  is 
known  from  spectroscopic  measurements  (3)  to  be  about  500  cal  ./mole. 

Equation  (7)  has  a  number  of  interesting  properties  which  are 
not  immediately  obvious.  When  Z  ■  Zy  *  h  and  x  *  5  or  less.  As  is 
greater  than  zero  for  all  positive  values  of  d/T.  Nevertheless, 
reasonable  values  of  d/T0  are  found  for  x  *  6  or  greater.  Positive 
values  of  T^  cannot  be  obtained  for  pentane  (x  •  5)  even  if  the  free 
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volume  Is  decreased  to  zero.  As  a  matter  of  fact,  adjustment  of  v 
does  little  to  Improve  the  fit  to  the  data  above  x  =  5,  and  It  seems 
expedient  to  eliminate  this  parameter  from  the  theory. 

It  was  noted  that  for  long  chains,  increasing  Z  decreases  the 
entropy.  Unfortunately,  for  x  »  5  the  trend  is  reversed,  as  shown 
in  Table  5,  where  the  terms  in  equation  (7)  involving  Z  are 
evaluated  for  v  ■  o.  Coordination  numbers  less  than  4  are  physically 
unreasonable,  so  the  best  choice  seems  to  be  4,  which  is  the  number 
used  by  Gibbs  and  DIMarzio. 


Table  5 


EFFECT  OF  COORDINATION  NUMBER  ON  CONFIGURATIONAL  ENTROPY 


z 

X 

Z-2  ,  (Z-2)  x  ♦  2 

2  Z  x 

x*1  In  /7z-2)  x  4  27  (Z-l)/2 

£S/R 

4 

5 

-0.511 

0.578 

0.067 

4 

00 

-0.693 

0 

-0.693 

6 

5 

-0.620 

0.801 

0.181 

6 

oo 

-0.81' 

0 

-0.811 

8 

5 

j69 

0.944 

0.275 

8 

OO 

.,.864 

0 

-0.864 

One  obvious  defect  in  the  Gibbs -DIMarzio  theory  is  the  failure  to 
take  into  account  the  fact  that  rotations  about  successive  bonds  are 
interdependent,  as  discussed  in  the  preceding  monthly  letter.  The 
mathematical  technique  for  handling  this  problem  was  developed  by 
LifsonC2^),  and  applied  to  the  calculation  of  random  coil  dimensions 
by  Nagal  and  Ishikawal2^)  and  Hoeve^2®),  and  to  the  entropy  of  melting 
oy  Starkweather  and  Boyd'2?).  The  decrease  in  entropy  due  to  this 
effect  is  at  a  maximum  for  £  /?  *  0,  and  becomes  negligible  for 
6/t'Z  5v25#27).  The  formula  given  by  Taylor'2®)  for  the  fraction  of 
"forbidden"  configurations  yields  0.05  R  for  the  decrease  in  entropy 
per  atom  in  pentane,  and  even  this  maximum  value  is  too  small  to  resolve 
the  difficulty  with  the  calculation  of  TQ .  For  the  infinite  polymer 
6/ T0^  3*9  (for  v  ■  o  and  Z  ■  Z/  *  4)  and  this  increases  as  x  decreases:, 
so  it  appears  that  modification  of  the  Gibbs-DiMarzio  theory  to  include 
the  pentane  effect  is  hardly  worthwhile. 

Starkweather  and  Boyd(2?)  point  out  that  lattice  treatments  are 
overly  restrictive  and  underestimate  the  entropy  of  liquids.  They 
propose  addition  of  a  contribution  due  to  "long  range  disorder"  on  the 
order  of  0.5  R.  However,  this  would  lead  to  negative  values  of  T0 
for  all  the  n-alkanes  below  elcosane. 
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For  the  reasons  given  above,  free  volume  and  long  range  disorder 
were  omitted,  and  Z  was  set  equal  to  Z/  «  4,  to  arrive  at  the  final 
expression  to  be  used  for  comparison  with  experiment: 


0=  *n  jHr  +  "x 

The  energy  difference  between  rotational  isomers  is  generally 
lower  in  the  liquid  state  than  in  the  gas,  hence  the  liquid  state 
value  should  be  used  in  equation  (8) .  Data  of  this  type  are  not 
plentiful,  so  the  most  practical  way  to  estimate  £  for  polymers  is 
by  application  of  equation  (8)  to  experimental  measurements  of  TQ. 

The  Mathematical  Analysis  Group  estimated  the  best  values  of  £ 
in  equation  (8)  for  the  set  of  T0(ty)  values  listed  in  Table  1 
(smoothed  data)  and  2  (raw  data)  for  kinematic  viscosities  of  the 
n -alkanes .  These  and  the  corresponding  TQ  calculated  for  the  infinite 
polymer  are  shown  in  Table  6.  It  was  pointed  out  above  that  positive 
values  of  £/Tq  are  not  obtained  for  x  ■  5  or  less,  so  these  data 
were  omitted  from  the  analysis.  Furthermore,  most  of  the  error  in 
estimation  was  contributed  by  hexane,  x  ■  6,  so  a  value  is  also 
shown  with  hexane  omitted. 


Table  6 


ROTATIONAL  ENERGY  TOR  n -ALKANES 


Smoothed  data 

Smoothed  data  minus  hexane 
Raw  data 


£  (cal  ./mole)  Tn(*°)(*K) 


493.4  129.5 
490.8  128.8 
467-4  122.6 


The  calculated  energy  differences  are  in  excellent  agreement 
with  the  spectroscopically  measured  value 5 (23)  0f  450  for  pentane, 
and  470  and  520  for  two  rotational  isomers  of  hexane.  Butane,  for 
some  reason,  has  a  higher  energy  difference,  76O  cai./mole,  than 
the  higher  alkanes . 

The  limiting  values  of  T0  for  x  ■  00  all  fall  within  a  1 
range,  and  this  seems  to  correspond  to  the  well-known  transition 
at  about  120*K'^^  in  polyethylene. 


In  view  of  the  success  of  the  Gibbs -DiMarzio  theory  in  pre¬ 
dicting  the  magnitude  of  T0,  it  becomes  of  pressing  importance 
to  find  a  correlation  of  6  with  the  structure  of  the  repeating 
unit.  The  difficulty  of  this  problem  is  illustrated  by  the  situ¬ 
ation  in  the  methacrylate  series,  where  the  Tg  of  poly(methyl 
methacrylate) ,  the  parent  member,  varies  over  a  range  of  at  least 
tiO*  depending  on  stereoregularity.  Simple  additive  schemes  are 
basically  incapable  of  accounting  for  6uch  effects,  and  detailed 
examination  of  the  molecular  structure  will  be  required. 

With  the  above  considerations  in  mind,  an  attempt  was  made 
to  devise  a  scheme  for  representing  diagranmntically  the  spatial 
requirements  of  polymeric  chains,  in  hopes  of  demonstrating  the 
steric  effects  which  influence  the  behavior  of  polymers  in  bulk. 

In  view  of  the  variations  in  bond  lengths  and  angles  in  different 
structures,  and  the  large  number  of  conformations  available  to 
each  monomer  unit,  it  is  necessary  to  use  some  simplifying 
assunptions  to  make  the  problem  tractable.  Furthermore,  this 
type  of  analysis  is  to  be  applied  to  polymers  which  have  not  yet 
been  synthesized,  so  that  knowledge  of  the  crystal  structure  cannot 
be  used  as  w  basis. 

The  most  satisfactory  procedure  found  so  far  is  as  follows: 
the  atoms  are  located  on  a  tetrahedral  lattice,  shown  in  pro¬ 
jection,  with  the  relative  levels  of  adjacent  atoms  depicted  by 
means  of  the  "flying  wedge"  symbol.  Bonds  between  superimposed 
atoms  therefore  appear  as  in  a  Newman  projection.  Superimposed 
atoms  are  shown  by  circles  of  different  sizes,  small  circles 
lying  above  the  larger  ones .  Figure  3  illustrates  the  use  of 
these  diagrams  in  evaluating  steric  effects  in  the  n-alkanes. 
Structures  I -IV  represent  pentane,  and  V  and  VI  are  sections  of 
a  longer  chain  of  methylene  units. 

The  use  of  a  tetrahedral  lattice  limits  the  conformation  of 
each  bond  to  the  three  isomers:  trans,  gauche  (left)  or  gauche 
(right),  labeled  T,  G  and  G7 .  Two  types  of  steric  interaction 
can  be  discerned.  Hydrogen  atoms  may  occupy  adjacent  lattice 
sites  (structures  II  and  III),  which  places  them  within  their 
van  der  Waais  radius  of  interaction,  o1"  may  occupy  the  same  site 
(structure  IV),  presumably  a  "forbidden"  configuration.  These 
atoms  are  l,o  to  each  other  in  the  first  case  and  1,7  in  the 
second,  arranged  in  the  form  of  a  cyclohexane  ring,  which  can 
appear  either  flat  or  on  edge.  Interactions  with  more  distant 
atoms  can  probably  be  neglected. 

As  Figure  3  shows,  the  all  trans  configuration  is  the  most 
stable,  while  the  sequence  G  G7  is  forbidden.  Since  forbidden 
configurations  occur  first  in  pentar.e,  Taylor named  it  the 
"pentane  effect"  . 
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Figure  3»  Projections  of  Molecules  on  Diamond  Lattice 


Substituents  on  the  backbone  can  lead  to  double  occupation  of 
lattice  sites  in  the  all  trans  configuration,  so  that  helical 
configurations  tend  to  be  favored.  Two  helical  chains  are  depicted 
in  Figure  3  to  illustrate  the  manner  in  which  they  appear  on  this 
type  of  diagram.  Structure  V  has  six  atoms  per  turn,  and  structure 
VI  has  four.  When  the  backbone  carries  substituents  it  is  not 
always  possible  to  find  a  configuration  free  of  the  pentane  effect. 
The  important  factor  in  chain  stiffness  is  presumably  the  occurrence 
of  one  configuration  which  is  strongly  preferred  over  all  others, 
although  it  may  contain  some  double  occupancies. 

It  will  be  necessary  to  examine  a  number  of  selected  polymers 
for  which  6.  has  been  determined  by  use  of  the  Gibbs -DiMarzio  theory, 
to  see  if  interactions  are  evident  on  the  diagrammatic  projections 
which  can  account  for  the  variations  in  stiffness  energy.  If  this 
approach  is  successful,  each  type  of  interaction  could  be  assigned 
an  energy,  depending  on  the  nature  of  the  atoms  involved,  and  a 
small  number  of  such  values  would  suffice  for  the  prediction  of  TQ 
for  a  large  number  of  different  polymers. 

Effect  of  Pressure  on  Viscosity 

Although  a  number  of  studies  have  been  made  of  the  effect  of 
pressure  on  viscosities  of  liquids,  as  yet  no  relation  has  been 
found  which  satisfactorily  represents  the  isotherms  observed.  It 
would  be  highly  desirable  to  be  able  to  separate  the  phenomenological 
effects  of  both  temperature  and  pressure  from  the  structural  and 
molecular  factors  involved  in  the  viscoelastic  behavior  of  polymers. 
This  in  fact  is  an  essential  first  step  in  the  correlation  of 
structure  with  properties.  An  approach  to  the  solution  of  this 
problem  is  described  below. 

The  Doolittle^)  equation  can  be  written 

In  9  ■  In  A^  ♦  bvg/vf  (9) 

where  Vf/vs  i8  the  relative  free  volume,  and  b  is  a  constant 
independent  of  temperature.  Comparison  with  equation  (1)  shows 
that  must  be  linear  in  T,  when  equation  (1)  is  obeyed.  Corres- 

V8 

pondingly,  the  influence  of  pressure  can  be  explained  in  terms  of 
a  reduction  in  free  volume  as  the  pressure  is  raised.  If  it  is 
assumed  that  the  "occupied"  volume,  v_,  i6  independent  of  pressure, 
then  a  suitable  equation  of  state  will  predict  the  effect  of  pressure 
on  free  volume  and  hence  viscosity.  Using  the  Tait  equation  and 
converting  equation  (9)  to  common  logarithms,  we  obtain 

log  ^  ■  log  A^  ♦  O.U343  bvg  fvf(l)  -  Kvx  log  (L  ♦  (10) 
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where  vj(l)  is  the  free  volume  when  no  pressure  is  applied,  is 
the  corresponding  specific  volume,  amd  K  and  L  are  the  constants 
in  the  Tait  equation.  Inverting  and  substituting  the  value  of 
vf(l)/bvs  from  equation  (1)  and  (9)#  **  have 


_ 

lo<j(rj/A^ 


T-T. 


2  3o\  Kv, 
bvs 


L  -  P 
L 


(11) 


K. is. known  to  be  independent  of  temperature,  and  in  fact  equals 
0.2058(^'  for  all  liquids  for  which  it  has  been  determined.  It  is 
evident  from  the  manner  in  which  equation  (11)  was  derived  that  A, 

B  and  TQ  are  independent  of  pressure.  The  change  in  viscosity  with 
pressure  is  therefore  predicted  by  means  of  only  two  adjustable 
parameters:  (bvg)  and  L.  In  accordance  with  the  finding  that 

equation  (1)  is  obeyed  more  closely  by  kinematic  than  absolute 
viscosities,  equation  (11)  has  been  applied  to  the  kinematic  vis¬ 
cosities  calculated  from  the  data  reported  by  Lovitz^l)  et  al.  It 
was  found  that  the  plots  of  l/log(n/A)  were  linear  in  log(L+P)/L 
for  every  liquid  examined,  for  the  proper  choice  of  L. 


The  free  volume  vf(l),  occupied  volune  vg  and  Doolittle 
constant  b  can  calculated  from  the  slope  and  intercept  of  the 
latter  plot.  These  are  shown  in  Table  7  along  with  the  constants 
A,  B  and  TQ  and  the  Tait  equation  parameter  L.  The  temperature- 
dependent  values  were  obtained  at  37.78",  except  for  n-octane  at 
30*. 


Table  7 


FREE  VOLUME 

PARAMETERS 

FOR  EFFECT  OF  PRESSURE 

ON  VISCOSITY 

Dodeca- 

hydrochrysene 

Liquid 

n-Octane 

9- ( 2-cyclohexylethyl ) 
hep tade cane 

-  Ferhydro- 
chrysene 

T0  (*K.) 

62.9 

155 

203 

250 

»,(**.) 

256.7 

370.7 

240.5 

215.7 

log  A, 

-1.634 

-1-3053 

-0.815 

-0.8616 

L  (atm) 

7U5 

700 

1200 

1200 

vf(l)  (cc/gr) 

0.5259 

0.4413 

0.1261 

0.1003 

vs  (cc/gr) 

0.9140 

0.7758 

0.9056 

0.8633 

vx  (cc/gr) 

1.4399 

1.2171 

1.0317 

0.9636 

b 

1-751 

3-229 

0.6180 

0-9530 

bvs  (cc/gr) 

i.600b 

2.505 

0.5597 

0.8227 

1)1  (cp) 

0.463 

12.1 

25.6 

514 

lb 

V 


There  are  seme  consistent  trends  4.n  these  results.  In  order 
of  increasing  T0:  Vf(l)  decreases,  v^  decreases  and  ^  increases. 
These  parameters  will  be  calculated  for  those  liquids  for  which 
viscosity-pressure  data  are  available.  It  is  hoped  that  the  data 
will  fit  into  either  a  corresponding  states  relationship  or  a  direct 
correlation  with  molecular  structure. 

Cutler,  et  al.^O)  reported  the  values  of  L  for  compressibilities 
of  some  of  their  compounds.  They  found  L  »  1275  atm.  for  9“(2- 
cyclohexylethyl)-heptadecane  at  37-78*,  whereas  the  viscosity  data 
fits  equation  (11)  only  for  L  *  700  atm.  The  use  of  absolute  instead 
of  kinematic  viscosity  in  equation  (11)  does  not  change  this  value 
of  L  appreciably.  The  L  values  for  viscosities  and  compressibilities 
will  be  compared  for  various  liquids,  to  see  if  this  situation  is 
general . 
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